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^ ' Abstract 



The excitation of soft dipole modes in light nuclei via inelastic electron scattering is investigated. 



I show that, under the proposed conditions of the forthcoming electron- ion colliders, the scattering 
cross sections have a direct relation to the scattering by real photons. The advantages of electron 
scattering over other electromagnetic probes is explored. The response functions for direct breakup 
are studied with few-body models. The dependence upon final state interactions is discussed. 
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A comparison between direct breakup and collective models is performed. The results of this 

in 

investigation are important for the planned electron-ion colliders at the GSI and RIKEN facilities. 
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I. INTRODUCTION 



Reactions with radioactive beams have attracted great experimental and theoretical in- 

I 

terest during the last two decades Progresses of this scientific adventure were reported 
on measurements of nuclear sizes [2J, the use of secondary radioactive beams to obtain in- 
formation on reactions of astrophysical interest 3, \4\, fusion reactions with neutron-rich 

n n it 

nuclei |5|, |6j, tests of fundamental interactions |7[, dependence of the equation of state of 



nuclear matter upon the asymmetry energy [8(, and many other research directions. Studies 
of the structure and stability of nuclei with extreme isospin values provide new insights into 
every aspect of the nuclear many-body problem. In neutron-rich nuclei far from the valley 
of /3-stability, in particular, new shell structures occur as a result of the modification of 
the effective nuclear potential. Neutron density distributions become very diffuse and the 
phenomenon of the evolution of the neutron skin and, in some cases, the neutron halo have 
been observed. 

New research areas with nuclei far from the stability line will become possible with 
newly proposed experimental facilities. Among these we quote the FAIR facility at the GSI 
laboratory in Germany. One of the projects for this new facility is the study of electron 
scattering off unstable nuclei in an electron-ion collider mode [9|. A similar proposal exists 



for the RIKEN laboratory facility in Japan 



10j. By means of elastic electron scattering, 



;hese facilities will become the main tool to probe the charge distribution of unstable nuclei 
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12j |. This will complement studies of matter distribution which have been performed 



in other radioactive beam facilities using hadronic probes. Inelastic electron scattering will 
test the nuclear response to electromagnetic fields. 

These facilities will provide accurate measurements of many nuclear properties of unstable 
nuclei. The reason is that electron scattering is a very clean probe. Its electromagnetic 
interaction with the nucleus is well understood. Inelastic electron scattering can also be 
very well described in the Born approximation. Higher order processes are only relevant for 
the distortion of the electron wavef unctions, affecting mostly electron scattering on heavy 
nuclei. 

Up to now, the electromagnetic response of unstable nuclei far from the stability line has 
)een studied with Coulomb excitation of radioactive beams impinging on a heavy target 
4)- This method has been very useful in determining the electromagnetic response in light 



2 



nuclei 



131 ]. For neutron- rich isotopes [14| the resulting photo- neutron cross sections are 



characterized by a pronounced concentration of low-lying El strength. The onset of low- 
lying El strength has been observed not only in exotic nuclei with a large neutron excess, 
but also in stable nuclei with moderate proton-neutron asymmetry. The problem with such 
experiments is that the probe is not very clean. It is well known that the nuclear interaction 
between projectile and target as well as the long range Coulomb distortion of the energy 
of the fragments interacting with the target (see, e.g. ref. are problems of a difficult 

nature. The nuclear response probed with electron does not suffer from these inconveniences. 

The interpretation of the low-lying El strength in neutron-rich nuclei engendered a de- 
bate: are these "soft dipole modes" just a manifestation of the loosely-bound character 



of 



IS 



ight neutron-rich nuclei, or are they a manifestation of the excitation of a resonance"' 
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19] . As far as I know, there has not been a definite answer to this simple ques- 



tion. This apparently innocuous question has nonetheless become the center of a even more 
widespread debate. It is believed that the weak binding of outermost neutrons gives rise 
to a direct break up of the nucleus and a consequent concentration of the electromagnetic 
response at low energies. The same weak binding can also lead to soft collective modes. 
In particular, the pygmy dipole resonance (PR), i.e. the resonant oscillation of the weakly- 
bound neutron mantle against the isospin saturated proton-neutron core. Its structure, 
however, remains very much under discussion. The electromagnetic response of light nuclei, 
leading to their dissociation, has a direct connection with the nuclear physics needed in sev- 



eral astrophysical sites [3j, |4j, [15|]. In fact, it has been shown [20( that the existence of pygmy 
resonances have important implications on theoretical predictions of radiative neutron cap- 
ture rates in the r-process nucleosynthesis, and consequently to the calculated elemental 
abundance distribution in the universe. 

In this work I study the general features of inelastic electron scattering off light nuclei, in 
particular their response in the continuum. An assessment of the theory of inelastic electron 
scattering appropriate for the conditions of electron-ion colliders is presented in section 
2. Special emphasis is put on the connection of electron scattering and the scattering by 
real photons, which will be useful to relate electron scattering and Coulomb dissociation 
measurements. Section 3 deals with the nuclear response within two and three-body models 
and their dependence upon final state interactions. Section 4 discusses the aspects of low 
energy collective modes in halo nuclei and their connection with the response obtained with 



3 



few-body models. The summary and conclusions will be presented in section 5. 



II. INELASTIC ELECTRON SCATTERING 



In the plane wave Born approximation (PWBA) the cross section for inelastic electron 



scattering is given by 
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+ 



Jhcf 
EE' - 



EE' + c 2 p • p' + m 2 c 4 



CL)f 



c 2 (p ■ q) (p' ■ q) — ^ 2 c 4 



F lj (q;ML)\ 2 + \F l ,(q;EL)\ 2 ] 



(1) 



c 2 (q 2 - qlf 

where Ji (Jf) is the initial (final) angular momentum of the nucleus, (E, p) and (£", p') are 
the initial and final energy and momentum of the electron, and (go, q) = ( ^ ^ > ^ P "^ P ^ ) 
is the energy and momentum transfer in the reaction. Fij (q; ILL) are form factors for mo- 
mentum transfer q and for Coulomb (C), electric (E) and magnetic (M) multipolarities, 
II = C, E, M, respectively. 

Here we will only treat electric multipole transitions. Moreover, we will treat low energy 
excitations such that E, E' ^> hajo, which is a good approximation for electron energies 
E ~ 500 MeV and small excitation energies AE = hcq® ~ 1 — 10 MeV. These are typical 



values involved in the dissociation o: 
Using the Siegert's theorem [23 



nuclei far from the stability line. 



241 ]. one can show that the Coulomb and electric form 



factors in eq. [T] are proportional to each other. Moreover, for very forward scattering angles 
(0 <C 1) the PWBA cross section, eq. [lj can be rewritten as 



da 



dVtdE^ 



^ dN { e EL) (E, E. 
\ dQdE^ 



7 , 



a 



(EL) 



(EL) 



7/ ' 



(2) 



where a\ (E^), with E y = hcqo, is the photo-nuclear cross section for the .EL-multipolarity, 
given by [4| 



r (EL) 
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a-: (E 1 



(27r) d (L + l) fE, 



2L-1 



dB (EL) 



(3) 



L[(2L + l)!!] 2 \hcj dE~ t 
In the long-wavelength approximation, the response function, dB (EL) /dE^, in eq. [3] is 
given by 



dB (EL) \(J f \\Y L (r)\\Ji. 



|2 r 



2X + 1 



dr r 



2+L 



Spf f L) (r) 



-i 2 



w(E. 



7/ ' 
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where w (E 7 ) is the density of final states (for nuclear excitations into the continuum) with 

(FT ) 

energy E 1 = Ef — Ei. The the transition density 5p\t (r) will depend upon the nuclear 
model adopted. 

For L > 1 one obtains from eq. (1) that 

~l 2L-1 



dNi EL) 



(E,E 



dVtdE^ 



x 



41 a 2E (6 
L + 1E [E^ Sm \2 

cos 2 (6/2) sin" 3 (0/2) 
1 + (2E/M A c 2 ) sin 2 (6/2) 
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L 



sm 



2 +tm 2 



One can also define a differential cross section integrated over angles. Since a\ does not 
depend on the scattering angle, this can be obtained from eq. Elby integrating dNe EL ^ / dttdE 1 
over angles, from 6 min = E^/E to a maximum value 6 m , which depends upon the experi- 
mental setup. 

Eqs. I2H5] show that under the conditions of the proposed electron-ion colliders, electron 
scattering will offer the same information as excitations induced by real photons. The 
reaction dynamics information is contained in the virtual photon spectrum of eq. El while the 
nuclear response dynamics information will be contained in eq. HI This is akin to a method 



(5) 



developed 
method 



ong time ago by Fermi [25j and usually known as the Weizsaecker- Williams 
26l |. The quantities dNe E / dQdE^ can be interpreted as the number of equivalent 
(real) photons incident on the nucleus per unit scattering angle Q and per unit photon energy 
Ej. Note that £"0 (monopole) transitions do not appear in this formalism. As immediately 
inferred from eq. HI for L = the response function dB (EL) / dE^ vanishes because the 
volume integral of the transition density also vanishes in the long-wavelength approximation. 
But for larger scattering angles the Coulomb multipole matrix elements (CL) in eq. [Hare in 
general larger than the electric (EL) multipoles, and monopole transitions become relevant 



27] . 



In figured] we show the virtual photon spectrum for the El, E2 and E3 multipolarities 
for electron scattering off arbitrary nuclei at E e = 100 MeV. These spectra have been 
obtained assuming a maximum scattering angle of 5 degrees. An evident feature deduced 
from this figure is that the spectrum increases rapidly with decreasing energies. Also, at 



excitation energies of 1 MeV, the spectrum yields the ratios dNe^/dN, 



(El) 
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FIG. 1: (Color online) Virtual photon spectrum for the El, E2 and E3 multipolar ities in electron 
scattering off arbitrary nuclei at E e = 100 MeV and maximum scattering angle of 5 degrees. 

dNe E3 ^ I 'dNe — 100. However, although dNe EL ^ /dE^ increases with the multipolarity 
L, the nuclear response decreases rapidly with L, and El excitations tend to dominate 
the reaction. For larger electron energies the ratios /N^ and /N^ decrease 
rapidly. 

Note that a similar relationship as eq. [2] also exists for Coulomb excitation Q] in heavy 
ion scattering. In figure [2] we show a comparison between the El virtual photon spectrum, 
dN e /dEy, of 1 GeV electrons with the spectrum generated by 1 GeV/nucleon heavy ion 
projectiles. In the case of Coulomb excitation, the virtual photon spectrum was calculated 
in ref. 4[, eq. 2.5.5a. For simplicity, we use for the strong interaction distance R — 10 
fm. The spectrum for the heavy ion case is much larger than that of the electron for large 
projectile charges. For 208 Pb projectiles it can be of the order of 1000 times larger than that 
of an electron of the same energy. As a natural consequence, reaction rates for Coulomb 
excitation are larger than for electron excitation. But electrons have the advantage of being 
a clean electromagnetic probe, while Coulomb excitation at high energies needs a detailed 
theoretical analysis of the data due to contamination by nuclear excitation. As one observes 
in figure [21 the virtual spectrum for the electron contains more hard photons, i.e. the 
spectrum decreases slower with photon energy than the heavy ion photon spectrum. This is 
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FIG. 2: (Color online) Comparison between the virtual photon spectrum of 1 GeV electrons (dashed 
line), and the spectrum generated by a 1 GeV/nucleon heavy ion projectile (solid line) for the El 
multipolar ity, as a function of the photon energy. The virtual photon spectrum for the ion has 
been divided by the square of its charge number. 

because, in both situations, the rate at which the spectrum decreases depends on the ratio 
of the projectile kinetic energy to its rest mass, E/mc 2 , which is much larger for the electron 
(m = m e ) than for the heavy ion (m = nuclear mass). 

To obtain an effective luminosity per unit energy, the equivalent photon number is mul- 
tiplied by the experimental luminosity, L e A, i.e. d ^ ff = L e A-^-- The number of events per 
unit time, N T , is given by the integral N T = f a(E 7 )dL e ff, where a-(E y ) is the photonu- 
clear cross section. Assuming that the photonuclear cross section peaks at energy E and 
using the Thomas- Reiche-Kuhn (TRK) sum rule we can approximate this integral by 
N T = dL d ^ x 6 x lCT 26 ^, where dL e ff/dE is expressed in units of cm~ 2 s~ 1 MeV~ 1 . The 
giant resonances exhaust most part of the TRK sum rule and occur in nuclei at energies 
around E = 15 MeV. For 1 GeV electrons dN(E 7 = E )/dE 7 « 6 x lfr 3 /MeV. With a 
luminosity of L e ^ = 10 25 cm _2 s _1 , one gets d ^ ff ~ 6 x 10 22 cm~ 2 MeV _1 s _1 and a number 
of events N T m 4 x 10~ 3 NZ/A w 10~ 3 A s _1 . Thus, for medium mass nuclei, one expects 
thousands of events per day. These estimates increase linearly with the accelerator luminos- 
ity, L e A_, and show that studies of giant resonances in neutron-rich nuclei is very promising 
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f the order of 5%-10%, of the TRK 



471 ] . However, these modes occur 



at the proposed facilities 0, Q. Only a small fraction, o 
sum-rule goes into the excitation of soft-dipole modes 
at a much lower energy, E r « 1 MeV, where the number of equivalent photons (see figure 
2) is at least one order of magnitude larger than for giant resonance energies. Therefore, 
inelastic processes leading to the excitation of soft dipole modes will be as abundant as 
those for excitation of giant resonances. However, one has to keep in mind that it is not 
clear if experiments with very short-lived nuclei will be feasible at the proposed electron-ion 
colliders. 



III. DISSOCIATION OF WEAKLY-BOUND SYSTEMS 



One-neutron halo 



In this section I will consider the dissociation of a weakly-bound (halo) nucleus from a 
bound state into a structureless continuum. I calculate the matrix elements for the response 



function in eq. H] with a two-body model which has 



29 



30 



3een used previously to study Coulomb 



31 



32 



33, 



34J. The initial wavefunction 



excitation of halo nuclei with relative success 

can be written as \& jm = r~ l uijj{r)yuM, where Rijj(r) is the radial wavefunction and 3^zjm 



is a spin-angle function 35|. The radial wavefunction, -u^j(r), can be obtained by solving 



the radial Schrodinger equation for a nuclear potential, Vjf- (r). Some analytical insight 
may be obtained using a simple Yukawa form for an s-wave initial wavefunction, Uq(t) = 
Aq exp(—i]r) , and a p-wave final wavefunction, u\(r) = ji(kr) cos<5i — rii(kr) sin#i. In 
these equations r\ is related to the neutron separation energy S n = h 2 r] 2 /2p, \i is the reduced 
mass of the neutron + core system, and hk = y/2pE r , with E r being the final energy of 
relative motion between the neutron and the core nucleus. A is the normalization constant 
of the initial wavefunction. The transition density is given by r 2 5pif (r) = e e ffAiUi(r)uf(r), 
where i and / indices include angular momentum dependence and e e ff = —eZ c /A is the 
effective charge of a neutron+core nucleus with charge Z c . The El transition integral 
Zy, = J °° dr r 3 5pif (r) for the wavefunctions described above yields 



1, 



2k 2 



s— >p 



r](rf + ?,k 2 y 
2P" 



(r] 2 + k 2 Y 
e e ffh 2 2E r 
2/x (S n + E r 



cos 5\ + sin 5\ 



\2h 2 ) -1/ax +pr 1 E r /h 2 _ 



(6) 
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FIG. 3: (Color online) |T s ^ p | 2 calculated using eq. [U assuming e e fj = e, A = 11 and S n = 0.5 
MeV, as a function of E r . The long dashed curve corresponds to a± = —10 fm -3 and r\ = —0.5 
fm _1 , the dashed curve corresponds to aq = —50 fm -3 and n = 1 fm -1 , the solid curve corresponds 
to a\ = r\ = 0, and finally, the dotted curve corresponds to oq = —10 fm -3 and r\ = 0.5 fm -1 . 

where the effective range expansion of the phase shift, k 2l+1 cot 5 ~ — l/a/ + r//c 2 /2, was used 
in the second line of the above equation. For / = 1, cq is the "scattering volume" (units 
of length 3 ) and rq is the "effective momentum" (units of 1/length). Their interpretation is 
not as simple as the / = effective range parameters. Typical values are, e.g. a\ = —13.82 
fm -3 and rq = —0.419 fm" 1 for n+ 4 He jq/2-wave scattering and cq = —62.95 fm -3 and 
r% = —0.882 fm" 1 for n+ 4 He p 3 / 2 -wave scattering [361 ] . 

The energy dependence of eq. |6] has some unique features. As shown in previous works 



29, 30 



371 ] , the matrix elements for electromagnetic response of weakly-bound nuclei present 
a small peak at low energies, due to the proximity of the bound state to the continuum. 
This peak is manifest in the response function of eq. |4] 

dB(EL) 2 £, L+1/2 m 

dE ° C| ^ 1 K (S n + E r ) 2L+2 - [ ) 

It appears centered at the energy j^O] Eq EL ^ ~ iq^/f ^n f° r a generic electric response of 
multipolarity L. For El excitations, the peak occurs at Eq ~ 3S n /5. 

The second term inside brackets in eq. [6] is a modification due to final state interactions. 
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This modification may become important, as shown in figure [3j where |T S ^ P | 2 calculated 
with eq. [6] is plotted as a function of E r . Here, for simplicity, I have assumed the values 
e e ff — e, A — 11 and S n = 0.5 MeV. This does not correspond to any known nucleus and is 
used to assess the effect of the scattering length and effective range in the transition matrix 
element. The long dashed curve corresponds to a\ = —10 fm -3 and r\ = —0.5 fm , the 
dashed curve corresponds to a\ = —50 fm -3 and r\ = 1 fm -1 , the solid curve corresponds 
to ai = ri = 0, and finally, the dotted curve corresponds to a x = —10 fm~ 3 and r\ = 0.5 
fm -1 . Although the effective range expansion is only valid for small values of E r , it is evident 
from the figure that the matrix element is very sensitive to the effective range expansion 
parameters. 

The strong dependence of the response function on the effective range expansion param- 
eters makes it an ideal tool to study the scattering properties of light nuclei which are of 
interest for nuclear astrophysics. It is important to notice that the one-halo has been stud- 
ied in many experiments, e.g. for the case of n Be for which there are many data available 



fsee refs. 



38 



39, 



401] ). In these papers one can find a detailed analysis of how the nuclear 



shell-model can explain the experimental data, by fitting the spectroscopic factors for several 
single-particle configurations. It is beyond the purpose of the present paper to reproduce 
theses data, in view of the simple model adopted above. The main goal of this section is to 
show the relevance of final state interactions. 



B. Two-neutron halo 

Many weakly-bound nuclei, like 6 He or 11 Li, require a three-body treatment in order to 
reproduce the electromagnetic response more accurately. In a popular three-body model, 
the bound-state wavefunction in the center of mass system is written as an expansion over 



hyperspherical harmonics (HH), see e.g. 411 ]. 



*(x,y) = ^ Yl ®K l Ls(p)\j l Ki(n 5 )®x S ] jM - (8) 



p 5/2 



x*>y 



Here x and y are Jacobi vectors where (see figure H|) x = 4g (r x — r 2 ) and y = 
^ 2(A-2) ( r i+ r a _ where A is the nuclear mass, ri and r 2 are the position of the nu- 
cleons, and r c is the position of the core. The hyperradius p determines the size of a 
three-body state: p 2 = x 2 + y 2 . The five angles {f^} include usual angles (0 x ,(f) x ), (9 y ,(py) 

10 




FIG. 4: (Color online) Jacobian coordinates (x and y) for a three-body system consisting of a core 
(c) and two nucleons (1 and 2). 

which parametrize the direction of the unit vectors x and y and the hyperangle 9, related 
by x = psin9 and y = pcos9, where < 9 < tt/2. 

The insertion of the three-body wavefunction, eq. El into the Schrodinger equation yields 
a set of coupled differential equations for the hyperradial wavefunction ^^ls 0°)- Assuming 
that the nuclear potentials between the three particles are known, this procedure yields the 
bound-state wavefunction for a three-body system with angular momentum J. 

In order to calculate the electric response we need the scattering wavefunctions in the 
three-body model to calculate the integrals in eq. HI One would have to use final wave- 



functions with given momenta, including their angular information. When t 
interaction is disregarded these wavefunctions are three-body plane waves 43, 



le final state 



44J. To carry 



out the calculations, the plane waves can be expanded in products of hyperspherical har- 
monics in coordinate and momentum spaces. However, since we are only interested in the 
energy dependence of the response function, we do not need directions of the momenta. 
Thus, instead of using plane waves, I will use a set of final states which just include the 
coordinate space and energy dependence. 



I will also adopt an approach closely related to the work of Pushkin et al. 



431 ] (see also 



45l . |46| ) . For weakly-bound systems having no bound subsystems the hyperradial functions 
entering the expansion [8] behave asymptotically as 42j $ a (p) — ► constant x exp (—rjp) as 
p — ► oo, where the two-nucleon separation energy is related to rj by S^n = h 2 rj 2 / (2m n). 
This wavefunction has similarities with the two-body case, when p is interpreted as the 
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FIG. 5: (Color online) Comparison between the calculation of the response function (in arbitrary 
units) with eqs. [9]and[10l using 5 nn = and 5 nc = 0, (dashed line), or including the effects of final 
state interactions (continuous line). The experimental data are from ref. [49(. 

distance r between the core and the two nucleons, treated as one single particle. But notice 
that the mass would have to be replaced by Im^ if a simple two-body (the dineutron- 



model 



371 ] ) were used for n Li or 6 He. 



Since only the core carries charge, in a three-body model the El transition operator is 



given by M ~ yY\u (y) for the final state (see also 44]]). The El transition matrix element 



is obtained by a sandwich of this operator between $ a (p) / p 5//2 and scattering wavef unctions. 



In ref. 



431 ] the scattering states were taken as plane waves. I will use distorted scattering 



states, leading to the expression 



1 (El) = ! dxdy ^ a ^ y 2 xu p (y) u q (x) , (9) 

where u p (y) = ji (py) cos 5 nc — ri\ (py) cos 5 nc is the core-neutron asymptotic continuum 
wavefunction, assumed to be a p-wave, and u q (x) = jo (qx) cos S nn — n (qx) cos S nn is the 
neutron-neutron asymptotic continuum wavefunction, assumed to be an s-wave. The relative 



momenta are given by q = ^ (qi — q 2 ), and p = y 2 ( A ~ 2 ^ ( kl + k2 — k c ) . 

The El strength function is proportional to the square of the matrix element in eq. 
[9] integrated over all momentum variables, except for the total continuum energy E r = 



12 



H 2 (q 2 + p 2 ) /2mN- This procedure gives 
dB (El) 



dE r 



constant x J \I (El)\ 2 E 2 cos 2 G sin 2 QdQdVL q dVt p , (10) 



where G = tan -1 (q/p). 

The 1 So phase shift in neutron-neutron scattering is remarkably well reproduced up to 
center of mass energy of order of 5 MeV by the first two terms in the effective-range expansion 
kcot5 nn ~ -1/a nn ~\~ r nnk 2 /2. Experimentally these parameters are determined to be a nn — 
—23.7 fm and r nn = 2.7 fm. The extremely large (negative) value of the scattering length 
implies that there is a virtual bound state in this channel very near zero energy. The 
■> wave scattering in the n- 9 Li ( 10 Li) system appears to have resonances at low energies 



481 ]. I assume that this phase-shift can be described by the resonance relationsin 5. 



(T/2)/yJ(E r - E R f + r 2 /4, with E R = 0.53 MeV and T = 0.5 MeV [48]. 

Most integrals in eqs. [9] and [10] can be done analytically, leaving two remaining integrals 
which can only be performed numerically. The result of the calculation is shown in figure [51 
The dashed line was obtained using 5 nn = and 5 nc = 0, that is, by neglecting final state 
interactions. The continuous curve includes the effects of final state interactions, with <L„ 



49|. The 



and 5 nc parametrized as described above. The experimental data are from ref. 
data and theoretical curves are given in arbitrary units. Although the experimental data is 
not perfectly described by either one of the results, it is clear that final state interactions 
are of extreme relevance. 



As pointed out in ref. 43|, the El three-body response function of n Li can still 



be described by an expression similar to eq. [7J but with different powers. Explicitly, 
dB (El) /dE r oc EH + £ r ) . Instead of S 2n , one has to use an effective S%~/ = aS2 n , 
with a ~ 1.5. With this approximation, the peak of the strength function in the three-body 
case is situated at about three times higher energy than for the two-body case, eq. [71 In 
the three-body model, the maximum is thus predicted at E 1 ^ 1 ^ ~ 1.852 n , which fits the 



experimentally determined peak position for the n Li El strength function very well 43]. 
It is thus apparent that the effect of three-body configurations is to widen and to shift the 
strength function dB (El) jdE to higher energies. 

It is worthwhile mentioning that the data presented in figure [5] and of other experiments 



16 



3 



al. 



is different in form and magnitude of the more recent experiment of Nakamura et 



51] . The reason for the discrepancy is attributed to an enhanced sensitivity in the 

13 



experiment of ref. [5JJ to low relative energies below E re i = 0.5 MeV compared to previous 
experiments. Also, this recent experiment agrees very well with the nn-correlated model of 
Esbensen and Bertsch [52|. This theoretical model is different than the model presented 
in this section in many aspects. In principle, the three-body models should be superior, as 
they include the interactions between the three-particles without any approximation. For 



example, ref. 52] use a simplified interaction between the two-neutrons. On the other hand, 



they include the many-body effects, e.g. the Pauli blockingof the occupied states in the 



core. It is not well known the reason why the data of ref. 5jj is better described with the 



model of ref. 52] than traditional 3-body models. 



IV. COLLECTIVE EXCITATIONS: THE PIGMY RESONANCE 
A. The hydro dynamical model 

We have seen that the energy position where the soft dipole response peaks depends 
upon the few body model adopted. Except for a two-body resonance in 10 Li, there was no 
reference to a resonance in the continuum. The peak in the response function can be simply 
explained by the fact that it has to grow from zero at low energies and return to zero at 
large energies. In few-body, or cluster, models, the form of the bound-state wavefunctions 
and the phase space in the continuum determine the position of the peak in the response 
function. Few-body resonances will lead to more peaks. 

Now I shall consider the case in which a collective resonance is present. As with giant 
dipole resonances (GDR) in stable nuclei, one believes that pygmy resonances at energies 
close to the threshold are present in halo, or neutron-rich, nuclei. This was proposed by 



Suzuki et al. 53( using the hydrodynamical model for collective vibrations. The possibility 
to explain the soft dipole modes (figure [5]) in terms of direct breakup, has made it very 
difficult to clearly identify the signature of pygmy resonances in light exotic nuclei. 

The hydrodynamical model, first suggested by Goldhaber and Teller 54] and by Stein- 
wedel and Jensen [55] needs adjustments to explain collective response in light, neutron- 
rich, nuclei. Because clusterization in light nuclei exists, not all neutrons and protons can 
be treated equally. The necessary modifications are straight-forward and discussed next. 
To my knowledge, the radial dependence of the transition densities in the hydrodynamical 
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FIG. 6: (Color online) Hydrodynamical model for collective nuclear vibrations in halo nuclei. The 
(a) Steinwedel-Jensen (SJ) mode and the (b) Goldhaber- Teller (GT) mode are shown separately. 

model for light, neutron-rich, nuclei has not been discussed in the literature. I will use the 
method of Myers et al. [56|], who considered collective vibrations in nuclei as an admixture 
of Goldhaber- Teller and Steinwedel-Jensen modes. 

When a collective vibration of protons against neutrons is present in a nucleus with charge 
(neutron) number Z (N), the neutron and proton fluids are displaced with respect to each 
other by d\ = a\R and each of the fluids are displaced from the origin (center of mass of 
the system) by d p = Nd\/A and d n = —Zdi/A. This leaves the center of mass fixed and 
one gets for the dipole moment D\ = Zed p = a\N ZeR/ A. The GT model assumes that the 
restoring force is due to the increase of the nuclear surface which leads to an extra energy 
proportional to A 2 / 3 . In this model, the inertia is proportional to A and the excitation 
energy is consequently given by E x oc a/ A 2 / 3 /A = A -1 / 6 . 

For light, weakly-bound nuclei, it is more appropriate to assume that the neutrons inside 
the core (A c , Z c ) vibrate in phase with the protons. The neutrons and protons in the core 
are tightly bound. An overall displacement among them requires energies of the order 
of 10-20 MeV, well above that of the soft dipole modes. The dipole moment becomes 
Di = edi(Z c A — ZA C )/A = Z^edi, where di is a vector connecting the center of mass of 
the two fluids (core and excess neutrons). We see that the dipole moment is now smaller 
than before because the effective charge changes from NZ/A in the case of the GDR to 
Z^)f = (Z C A — ZA C ) /A. This effective charge is zero if Z C A = ZA C and no pigmy resonance 
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FIG. 7: (Color online) Contour plot of the nuclear transition density in the hydrodynamical model 
consisting of a mixture of GT and SJ vibrations. The darker areas represent the larger values of 
the transition density in a nucleus which has an average radius represented by the dashed circle. 
The legend on the right displays the values of the transition density within each contour limit. 

is possible in this model, only the usual GDR. 

Figure [6] shows a schematic representation of the hydrodynamical model for collective 
nuclear vibrations in a halo nucleus, as considered here. Part (a) of the figure shows the 
Steinwedel- Jensen (S J) mode in which the total matter density of both the core and the halo 
nucleons do not change locally. Only the local ratio of the neutrons and protons changes. 
Part (b) of the figure shows a particular case of the Goldhaber- Teller (GT) mode, in which 
the core as a whole moves with respect to the halo nucleons. 

For spherically symmetric densities, the transition density in the GT mode can be calcu- 
lated from 5p p = p p (|r — d p |) — p p (r), where p p is the charge density. Using d\ <C R, it is 
straight-forward to show that 5p p ^ (r) = Spjp (r) Y w (r) , where 



and po is the ground state matter density. 

In the Steinwedel- Jensen (SJ) mode, the local variation of the density of protons is found 




(11) 
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(2) (21 

to be 5pp (r) = 5p p (r) Y w (r), where 

H 2) (0 = \l ^-Zf] } a 2 Kj x (kr) p (r) , (12) 



where K = 9.93. If the proton and neutron content of the core does not change 53j, the 
effective charge number in the SJ mode is given by Z^X = Z 2 (N — N C )/A(Z + N c ). 

The transition density at a point r from the center-of-mass of the nucleus is a combination 
of the SJ and GT distributions and is given by 5p p (r) = 5p p (r) Y w (r), where 

% (r) = \Iy R { Z %^Tt + Z ^f a2 Jl 31 ^}p^- ( 13 ) 
Changes can be accommodated in these expressions to account for the different radii of the 
proton and neutron densities. 

Figure [7] shows the contour plot, in arbitrary units, of the nuclear transition density in 
the hydrodynamical model, consisting of a mixture of GT and SJ vibrations. The darker 
areas represent the larger values of the transition density in a nucleus which has an average 
radius represented by the dashed circle. In this particular case, I have used the HF density 
jl2 . IotI for 11 Li, and a radius R — 3.1 fm. The parameters a\ and a 2 were chosen so that 
^e// a i = ^e// a 2, i-e. a symmetric mixture of the SJ and GT modes. 

Figure M shows the transition densities for n Li for three different assumptions of the 
SJ+GT admixtures, according to eq. [131 The dashed curve is for a GT oscillation mode, 
with the core vibrating against the halo neutrons, with effective charge number Z^X = 6/11, 
radius R = 3.1 fm, and a.\ — 1. The dotted curve is for an SJ oscillation mode, with effective 
charge number Z^X = 2/11, and a 2 = 1. The solid curve is their sum. Notice that the 
transition densities are peaked at the surface, but at a radius smaller than the adopted 
"rms" radius _R = 3.1 fm. 

The liquid drop model predicts an equal admixture of SJ+GT oscillation modes for large 
nuclei 56]. The contribution of the SJ oscillation mode decreases with decreasing mass 
number, i.e. a 2 — > as A — > 0. This is even more probable in the case of halo nuclei, 
where a special type of GT mode (oscillations of the core against the halo nucleons) is likely 
to be dominant. For this special collective motion an approach different than those used in 



refs. [56J and 
and Teller 



)3J has to be considered. The resonance energy formula derived by Goldhaber 
54j changes to 

EpR = [2aRm N A r ) ' (M) 
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FIG. 8: (Color online) Hydrodynamical transition densities for n Li and three different assumptions 
for the SJ+GT admixtures, according to eq. [13] The dashed curve is for a GT oscillation mode, 
with the core vibrating against the halo neutrons, with effective charge number -Z^// = 6/11, radius 
R = 3.1 fm, and a\ = 1. The dotted curve is for an SJ oscillation mode, with effective charge 

(2) 

number Z e L = 2/11, and a 2 = 1. The solid curve is their sum. 

where A r = A c (A — A c ) /A and a is the length within which the interaction between a 
neutron and a nucleus changes from a zero-value outside the nucleus to a high value inside, 
i.e. a is the size of the nuclear surface, ip is the energy needed to extract one neutron from 
the proton environment. 
Goldhaber and Teller 



54J ] argued that in a heavy stable nucleus if is not the binding energy 
of the nucleus, but the part of the potential energy due to the neutron-proton interaction. 
It is proportional to the asymmetry energy. In the case of weakly-bound nuclei this picture 
changes and it is more reasonable to associate (p to the separation energy of the valence 
neutrons, S. I will use if = /3S, with a parameter (3 which is expected to be of order of one. 
Since for halo nuclei the product aR is proportional to S" 1 , we obtain the proportionality 
E PR oc S. Using eq. [TUfor n Li , with a = 1 fm, R = 3 fm and <p = S 2n = 0.3 MeV, we 
get Epr = 1.3 MeV. Considering that the pygmy resonance will most probably decay by 
particle emission, one gets E r ~ 1 MeV for the kinetic energy of the fragments, which is 
within the right ballpark (see figure [5]). 

Both the direct dissociation model and the hydrodynamical model yield a bump in the 
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response function proportional to S, the valence nucleon(s) separation energy. In the direct 
dissociation model the width of the response function obviously depends on the separation 
energy. But it also depends on the nature of the model, i.e. if it is a two-body model, like 
the model often adopted for n Be or 8 B, or a three-body model, appropriate for n Li and 
6 He. In the two-body model the phase-space depends on energy as p(E) oc d 3 p/dE oc VE, 
while in the three-body model p (E) ~ E 2 . This explains why the peak of figure [5] is pushed 
toward higher energy values, as compared to the prediction of eq. It also explains the 
larger width of dB/dE obtained in three-body models. In the case of the pigmy resonance 
model, this question is completely open. 



The hydrodynamical model predicts [56| for the width of the collective mode V = fw/R, 
where v is the average velocity of the nucleons inside the nucleus. This relation can be 
derived by assuming that the collective vibration is damped by the incoherent collisions of 
the nucleons with the walls of the nuclear potential well during the vibration cycles (piston 



model). Using v = 3vp/4, where vp = \j2Ep Jtun is the Fermi velocity, with Ep = 35 MeV 
and R = 6 fm, one gets r ~ 6 MeV. This is the typical energy width a giant dipole resonance 
state in a heavy nucleus. In the case of neutron-rich light nuclei v is not well defined. There 
are two average velocities: one for the nucleons in the core, v c , and another for the nucleons 
in the skin, or halo, of the nuc 



\/v c v h- Following ref. 



eus, v/j. One is thus tempted to use a substitution in the form 
the width of momentum distributions of core fragments in 
knockout reactions, <r c , is related to the Fermi velocity of halo nucleons by yp = a/ 5cr 2 /m^. 
Using this expression with a c ~ 20 MeV/c, we get T = 5 MeV (with R = 3 fm). This value 
is also not in discordance with experiments (see figure [5]). 

Better microscopic models, e.g. those based on random phase approximation (RPA) 



calculations 



59 



60J are necessary to study pigmy resonances. The halo nucleons have 



to be treated in an special way to get the response at the right energy position, and with 



approximately the right width 571, [60] . Electron scattering will provide a unique opportunity 



to clarify this issue due to its better resolution over Coulomb excitation. 
B. Total inelastic cross sections 

One might argue that the total breakup cross section would be a good signature for dis- 
cerning direct dissociation versus the dissociation through the excitation of a pigmy collective 
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vibration. The trouble is that the energy weighted sum rule for both cases are approximately 



of the same magnitude 
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This can be shown by using the electric dipole strength 



function in the direct breakup model of ref. [37J, namely 

dB(El) _ n Zhe 2 Zl ff ^S- n (E - S n f 2 

where E = E r + S n is the total excitation energy. C is a constant of the order of one, 
accounting for the corrections to the wavefunction used in ref. 

The sum rule for dipole excitations, Si (El) = dE E dB ^ , is 

with Z 2 ff = (Z C A — ZA C ) 2 /[AA C (A — A c )]. This is the same (with C — 1) as eq. 1 of ref. 



471 ] , which is often quoted as the standard value to which models for the nuclear response in 



the region of pigmy resonance should be compared to. The response function in eq. [T5], with 
C = 1, therefore exhausts 100% of the so-called cluster sum rule 47J. The total cross section 
for electron breakup of weakly-bound systems is roughly proportional to Si. This assertion 
can be easily verified by using eqs. [2] and [151 assuming that the logarithmic dependence of 
the virtual photon numbers on the energy E = E 1 can be factored out of the integral in eq. 
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The dipole strength of the pigmy dipole resonance is given by the same equation [T6l 
The constant C is still of the order of unity, but not necessarily the same as in eq. [15] and 
the effective charges are also different. For the Goldhaber- Teller pigmy dipole model the 
effective charge is given by Z^h = [Z C A — ZA C ) /A, whereas for the Steinwedel- Jensen it is 
Z^jj = (Z 2 /A)(N — N c )/ (Z + N c ). Assuming that the Goldhaber- Teller mode prevails, one 
gets the simple prediction for the ratio between the cross sections for direct breakup versus 
excitation of a pigmy collective mode: 

^.direct A 

° = c - (17) 

a pigmy A c (A - A c ) V ' 

For n Li this ratio is 11C/18 while for n Be it is 11C/10. 



V. SUMMARY AND CONCLUSIONS 



I have studied the feasibility to determine low energy excitation properties of light, exotic, 
nuclei from experimental data on inelastic electron scattering. It was shown that for the 
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conditions attained in the electron-ion collider mode, the electron scattering cross sections 
are directly proportional to photonuclear processes with real photons. This proportionality 
is lost when larger scattering angles, and larger ratio of the excitation energy to the electron 
energy, E 1 /E 1 are involved. 

One of the important issues to be studied in future electron-ion colliders is the nuclear 
response at low energies. This response can be modeled in two ways: by a (a) direct breakup 
and by a (a) collective excitation. We have shown that in the case of direct breakup the 
response function will depend quite strongly on the final state interaction. This may become 
a very useful technique to obtain phase-shifts, or effective-range expansion parameters, of 
fragments far from the stability line. In the case of collective excitations, a variant of the 
Goldhaber- Teller and Steinwedel- Jensen model was used to obtain the transition densities 
in halo nuclei. 

The pygmy resonance lies above the neutron emission threshold, effectively precluding its 
observation in (7,7') experiments on very neutron-rich nuclei. Nonetheless, electron scatter- 
ing experiments will probe the response function under several conditions, including different 
bombarding energies, different scattering angles, etc. The study of pygmy resonances and of 
final state interactions will certainly be an important line of investigation in these facilities. 
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